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Abstract 

We consider a non-minimally coupled inflaton, in a higher order curvature background, leading 
to a potential which evolves with the curvature scalar of the Universe, and which describes 
two regimes. The first one is a de Sitter phase, where the potential is static, and an exact 
exponential solution is found for the inflaton. The second regime, triggered by the inflaton 
reaching a threshold, leads to a power-law expansion, during which the potential becomes flat, 
quickly enough for the inflaton never to reach the minimum of the initial symmetry breaking 
potential. This scenario is an alternative to the inflaton oscillating about a minimum of the 
potential, and where preheating is a consequence of the flattening of the potential during the 
lO ' power-law expanding phase. 

1 Introduction 

Among the numerous inflationary models, those of type II involve of a spontaneous symmetry 
breaking potential, with an inflaton rolling down the concave part of the potential, to end up 
in one of its a minima pQ. It is known that such concave potentials generate the so-called 
spinodal instability, where the absence of restoration force implies that fluctuations are not 
perturbative. As long as one neglects perturbations about an equilibrium state, the spinodal 
instability can be disregarded. But if one wishes to take into account quantum effects, such a 
concave potential is not realistic: it is known that the effective potential of a scalar theory is 
q ■ necessarily convex [T2] . 

It was suggested in [2] and in [3] that quantum fluctuations can be sufficient to restore 
symmetry, where this idea was also supported by numerical studies. In these works, spinodal 
instability is related to preheating, generated by the growth of quantum fluctuations. We will 
recover these ideas here, explaining that the restoration of symmetry is achieved by a flat 
■ effective potential. 

An explicit flattening of the potential is shown in [3], where the authors consider a descrip- 
tion based on the spinodal decomposition, using a potential function of two variables: the mean 
inflaton field and its fluctuations condensate, in the framework of a Hartree approximation. In 
this way, they obtain a potential which flattens as the condensate increases, taking thus into 
account the effects of quantum fluctuations. This work was also used in [5] for the study of 
particle production. 

We propose here another effective description, taking into account the time-dependence 
of the inflaton potential, where the latter evolves with the curvature scalar of the Universe, 
mimicking a Wilsonian coarse grained potential which evolves with the decreasing energy scale. 
Indeed, as was shown in [6], exact renormalization group transformations in the presence of a 
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spinodal instability lead to a flat effective potential which interpolates the minima of the bare 
potential. The resulting model is a non-minimally coupled inflaton, in a higher order curvature 
metric, similarly to other works already done |7J. The scenario we obtain contains then two 
phases: 

• Inflation is provided by a de Sitter phase, where we find an exact solution for the inflaton, 
which increases exponentially in time as long as it is in the spinodal region. During this 
inflationary period, the potential seen by the inflaton is static, since the curvature scalar 
is constant. We check that the backreaction of possible inflaton fluctuations on the metric 
are negligible; 

• The end of inflation occurs when the inflaton reaches a certain threshold to be defined. 
Assuming a power-law expansion, one finds that the scalar curvature decreases quickly 
enough for the potential to flatten before the inflaton reaches a minimum of the bare 
potential. As a consequence, the naive oscillations around the minimum do not occur, as 
was initially proposed in [8], and we show that this stage can indeed be identified with 
a preheating period, where the non-perturbative feature corresponds to the potential 
flattening in the spinodal region. 

The next section starts with a model of non-minimally coupled inflaton, with higher-order 
curvature terms. The non-minimal coupling of the inflaton with the metric plays the role of a 
potential for the inflaton. The former is therefore not static in general, since its characteristics 
depend on the metric. The Appendix explains some properties of the effective potential, in 
a Quantum Field Theory, and shows how the present potential was motivated, although the 
present study is based on classical Physics. 



2 Inflaton dynamics 

Motivated by the construction of a potential which can evolve with the geometric charac- 
teristics of the Universe, we assume the following symmetry breaking potential, based on a 
non-minimally coupled inflaton, 

if0 2 <0 2 p £7(0) = M( i2 _^)(2m a -ei2 + eA* !, )-|(«-A*V W 

if 2 > <P% U{<t>) = - imV + 1/ 

6 

with <£, = -(m 2 -£# + £^2), 
9 

where R is the curvature scalar of the Universe, and /U 2 ,£ are parameters to be determined 
later. The field <ft sp defines the limits of what will be called the spinodal region <fi 2 < <p 2 sp . The 
limit (f) sp is smaller than the minimum of the symmetry breaking potential, and varies with 
R. In the spinodal region, the non-minimal coupling plays the role of a quadratic potential 
for the inflaton, with a negative slope. The higher order curvature term is present in order to 
have a vanishing minimum for the potential. This potential is motivated by results in exact 
renormalization studies, as explained in the Appendix, and is sketched in figiU It leads to the 
following actions to study (m p i is the Plank mass): 
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• Inside the spinodal region, for (f> 2 < cj) 2 p : 

Sin = J d4 ^{(^-j( m2 + ^ 2 ) + l^) R 

- + f^ 2 - ^ 2 <P 2 + + ^ 2 )}, (2) 

such that the inflaton sees a potential which can depend on time, since this potential 
varies with R; 

• Outside the spinodal region, for (f> 2 > (f) 2 sp : 

Sout = I d'xV^^R - \ 9l jr4*r4> ~ + ~™V - ^> 4 } , (3) 

where the inflaton sees a static potential. 




2.1 Equations of motion 

From now on, we consider the metric ds 2 = ■ 
time only. The scalar curvature is then 

R = 6 



dt 2 + a 2 (t)dr 2 , and an inflaton which depends on 
'a d 2 \ 
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and a dot represents a time derivative. The actions to study are: 
• in the spinodal region 



Si n [a,4>] = J d 4 x a 3 |^(i?- \i 



\l67r g J 



which leads to the following equations of motion 

SS. 
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= <t> + 3H<t>-£(R-fj, 2 )<l>, (6) 



and 



' )S '" =^(2m 2 + ^ 2 )+2(^-^(m 2 + ^ 2 ))(2H + 3H 2 ) (7) 



2 \ J.2 



5a 2g " \1Qtt g 

+ ^0 2 + 2£ (0 2 + U + 2^00) + £ feti + 3# 2 - ^ 

+ — f 2# {3) + 12M + 9H 2 + 18H 2 H 
9 

where is the third time derivative of the Hubble parameter H = a/a] 
outside the spinodal region 

S out [aA} = Jdxa (^+-0 + - -0 j 



$S ou t n '1 , nrrl 2 1 , 9 i'A 



leading to the equations of motion: 

= + 3#0 - m 2 + ^0 a , (9) 

UlfJ o 

and 

^ = = ^ (2H + 3H 2 ) + -0 2 - — + im 2 2 - (1Q) 
5a 8tt V y 2 r 2# 2 ^ 24^ v ' 

In what follows, we will consider R > fi 2 , in order not to have a flat potential from the beginning. 
2.2 De Sitter inflation 

We look here for a solution where H = H = constant, such that R = R$ > [i 2 implies 
12H 2 > /i 2 . Eq.© gives 

+ 3tf o = £(12tf o 2 -/i 2 )0 (11) 
for which solutions are d> = $e Wl *, where 



ut = — ^- -1 ± 



V + f( 12 -w,n, (12) 



and uj^ > 0, whereas < 0. We are interested in the increasing solution, and therefore will 
disregard . 

To solve eq. (JTJ), we note that the first line of this equation is a contant, that we impose to 
vanish. This leads to the following constraint on the parameters of the model: 



3^ 

2g 



(2m 2 + £/i 2 ) + 



~8^° 



Eq-© gives then 



1 ■ 



2 + 2£ U 2 + 



for which solutions are 



$e^ 2 *, where 



^(m 2 + ^ 2 )tf 2 . 



< (3tf 2 --/i 2 



(13) 



(14) 



LOT, 



-1 ± 



\ 




(15) 



We can see that, if fi 2 > 6H^, 0J2 is positive and can be identified to uf. Together with the 
bound Rq > /i 2 , we impose then 



6H 2 < if < 12H { 



! 



and Co>2~ 



ujJ gives 



3(1 + 80 -1 + 
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i+_8_e ( 

16£ 



H 2 



(16) 



(17) 



which is to be understood as a condition the parameters £, fi 2 , H Q have to satisfy. For example, 
if /j, 2 = IOHq, we find numerically £ ~ 0.6. The condition f fl3|) leads then to 



10.8 Hi - 7.2 m 2 ~ ^-m 2 pl , 



;is) 



which is a reasonable assumption. 

We therefore showed here that an exponentially growing inflaton is an exact solution of the 
equations of motion, if we assume a de Sitter inflation, and no slow roll approximation was 
necessary. During this period of inflation, the potential seen by the inflaton is static, since 
R = Rq is constant. This regime is valid until the inflaton arrives at the border of the spinodal 



region, defined by 



sp 



6(m — £i?o + £/•* )l 9-, after which it sees the quartic potential 



outside the spinodal region, given in eqs.^TJ. The time to when inflation stops is then given by 



2„2u;+ 1 



6 2 
— [m 
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12£# 2 + £/i 2 ), 



(19) 



where we note u 



00 1 = u 2- The number of e-foldings the Universe has undertaken is then 

^ 2 12£tf 2 + £/, 2 \ 



N = H t 



Ho I ur 
m 

2u+ 



#$ 2 /6 



(20) 
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and can take many values, depending on the parameters of the model, but mainly on the initial 
amplitude If we consider the previous example, where /i 2 = IOHq and thus £ ~ 0.6, we 
obtain 

TV B 1.4 x In ( m ' -■ \f H ° V (21) 

such that we need m 2 > 1.2Hq, which is consistent with the assumption f fl8|) . since the latter 
implies m 2 < 1.5Hq. To fix the idea, we will take later m 2 ~ 1.3//q. 



2.3 Stability of the de Sitter phase 

The previous subsection is based on classical equations of motion for the inflaton and the 
metric, and we study here the stability of the de Sitter metric against possible fluctuations of 
the inflaton. Indeed, gravitational backreaction on a classical de Sitter background can have 
large effects (see [9], [10], and references therein, where conditions are derived for backreactions 
to be perturbative, and where studies beyond the first order of fluctuations are done). 
We start from the equation of motion fill I) , and consider its linear perturbation arising from 
the small fluctuations 8(f) and SH at a given time. Taking into account the solution = $e w ', 
where uj = oj± = o;^", such that 



5<p = u 2 S(j), 



(22) 



we obtain then 



SHn 



ai2H 2 - /i 2 



to 



3uiHn 



H 3H (u-8^H ) (23) 

Note that the second equation of motion fTl4|) does not give any more information, since it is 
quadratic in <p and, after differentiation and using eqs.f l22p . a division by <ft leads to the relation 
( I23|) . provided £ ~ 0.6, where the two equations of motion are consistent. If, for this value of 
£, we take /i 2 = IOHq, we obtain then uj ~ 0.357ifo, and the perturbations are related by 



SH 




5<p 


H 


= ei 






where e\ ~ 10 



(24) 



The action ((2]) is quadratic in the inflaton, such that one can assume that the spectral modes 
of 5<fi are uncorrelated and lead to the diffusion law \S(f)\ ~ H ^H t pi)]. We obtain then, from 
eq.(j2T 



6H n 



■wt 



(25) 



and if we use the number of e- foldings (1211) . with m 2 = 1.3Hq as stated above, this ratio is, at 
the end of inflation, 



5H n 



e 2 \l gN exp(e 3 N) } where e 2 ~ 1.3 x ei and e 3 ~ 1.4 x 10 , 



(26) 



end 



which is very small for a typical value like N = 70. It is therefore consistent to assume that 
backreactions of inflaton fluctuations on the metric can be neglected in this model. Note that 
this result is valid because: 
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• The non-minimal coupling in the action (J2J) gives rise to the term £(12£/q — fi 2 ) in the 
numerator of eq. (l23[) . which almost compensates the derivatives terms u 2 + 3u>Hq, and 
leads to the small parameter e\ and therefore €2', 

• The inflaton grows exponentially, and, at the end of inflation, the small factor e _tJi ° ~ 
e -o.357Af j n eq. f )25j) almost compensates the large factor H /Q oc \J ge N / lA , which leads to 
the small parameter €3. 

These two reasons are actually not independent, since the specific exponential solution for the 
inflaton is a result of the non-minimal coupling term in the original action. 



2.4 Preheating 

As stated before, inflation ends when the dilaton <fr reaches the limit of the spinodal region (f) sp , 
since the dilaton sees then another potential, given in action ([3]). The potential is continuous at 
the junction between the spinodal region and the outside, but differentiable once only, such that 
its second derivative, related to the effective mass of the dilaton, is discontinuous at <fi = <p sp . 
The transition between inflation and preheating is a consequence of this non-differentiable 
phenomenological potential. 

After exiting the spinodal region, the inflaton naively rolls down the minimum of the bare 
potential and oscillates about the minimum. This picture would be correct if the potential 
was static, but if the Universe enters a power-law expanding regime, the curvature scalar R 
decreases and the potential flatens, in the spinodal region. We will see that this flattening occurs 
much quicker than the inflaton needs to reach the minimum of the bare potential, such that 
the inflaton actually does not oscillates. A similar situation happens in the preheating scenario 
of [8], and we will see that this non-perturbative potential flattening can indeed correspond to 
a preheating period. 

To start with, let us review the naive oscillation mechanism. We assume small oscillations 
about the minimum of the bare potential, occurring at min given by <^ min = 6m 2 /g. The 
angular frequency of the oscillations is y2m, and we look for a solution of the form 

(f> ~ (f) min + At q cos(v / 2mt), (27) 

where At q << m j n . Assuming a power-law expanding Universe, with H = p/t, it is easy to 
see that the equations of motion ([2]) and (flOl) are satisfied to the first order in At Q /(p m i n if 

2 

p = - and q = —1. (28) 
3 

Let us come back now to the evolution of the potential inside the spinodal region. With the 
solution p = 2/3, the curvature scalar is 



R = — 2 i ( 29 ) 



4 
3f 

and decreases from Rq = 4/(3tg) until it reaches the value /1 2 . This takes the time 



Af ^G-vW' (30) 
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that we have to compare with the period of oscillations t = 27r/(V2m) of the solution (dTJ). For 
this, we come back to the example where /i 2 = IOHq, an d we consider the case m 2 = 1.3Hq, as 
suggested earlier. We find then 

At ~ O.Olr, (31) 

which means that the inflaton does not have time to oscillate: almost as soon as it exits the 
spinodal region, the potential flatens, which corresponds to a non-perturbative mechanism. 
We can now estimate the energy lost by the inflaton during this phase. At the end of inflation, 
the inflaton is at the border of the spinodal instability <fi = <p sp , such that the potential is 

U( ( j> sp ) = ^(R -fi 2 ) 2 , (32) 

and it vanishes at the time t + At. With our previous example fi 2 = IOHq, we find that the 
potential energy loss is then 

U(<f> sp ) ~ ytf 4 , (33) 

and, assuming an instant thermalization, the corresponding temperature is then 

T ~ g- 1/A H 0) (34) 

which can be sufficient to reach a GUT temperature ~ 10 _3 m p ;, if we assume the expected 
orders of magnitude g ~ 10 -12 and H ~ 10~ 6 m p ;. 



3 Concluding remarks 

As the potential inside the spinodal region becomes flat, quantum fluctuations of the inflaton 
become large, and can take any value between the minima of the initial potential. But the slope 
of the potential outside the spinodal region generates restoring forces, imposing oscillations of 
the inflaton inside the spinodal region. These large oscillations lead to the decay of the inflaton. 
As a consequence, the Universe is filled by the decay products resulting from this preheating 
period, and equations of motion involving the inflaton only are not valid anymore. What is 
happening after depends on these particles produced, and on their coupling to gravity, but this 
is left for a future work. 

We note that the spinodal instability, leading to a convex potential, is universal and indepen- 
dent of the precise shape of the potential. The present study was done with a quartic potential, 
but a cosine, as in the natural inflation scenario would lead to a similar mechanism. 

In the present work, backreaction of particles produced has been neglected, but taking these 
into account would not change the inflationary period described here, as well as the exiting of 
the spinodal region. The preheating period, though, might see its metric modified by the decay 
products, and this is a possible extension of this work. 

Finally, we note that, when R reaches /i 2 , the potential is flat and the Universe is filled 
with products of the inflaton decay, which, unlike the previous inflaton, are expected to be 
minimally coupled to gravity, and another phase of the Universe begins. 
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Appendix 



We review here some properties of the effective potential in a quantum field theory, leading 
to what motivated the model ([I]). We first show that the dressed potential is necessarily 
convex. Then we show the equivalence of the effective potential defined by the Legendre 
transform, which is the derivative-independent part of the proper graph generator functional, 
and the effective potential defined in the Wilsonian sens. Finally, we show how the coarse 
grained potential, defined by Wilsonian renormalization, interpolates the bare "Mexican-hat" 
potential with the dressed flat potential. This interpolation is obtained when the energy scale 
decreases, and the flat potential is obtained in the IR limit of the Wilsonian renormalization 
group transformations. 



3.1 Convexity 

We review here the construction of the effective potential, which contains the interactions 
dressed by quantum fluctuations, and we show that it is necessarily convex (12] . 
Since we are interested in the effective potential rather than the full effective action, we will 
consider a constant classical field and therefore define the partition function for a constant 
source j, which is, for a Euclidean metric, 

Z(j) = J V[4>] exp (-S[$\ - Vj$) , (35) 

where V is the volume of space time, and $ is the constant mode of the dynamical variable <fi. 
The connected graphs generator functional is W(j) = — In Z(j) and has derivative 



dW V 



dj 



- J V[$\ I exp (-S[4>] - Vj$) = V (!>) = V$, (36) 



where $ defines the corresponding constant classical field. It is easy to see that the second 
derivative of W is 

W m = ^ = v 2 ((<J) 2 - ($ 2 )) <0, (37) 

and is negative, for any value of the source j. The effective potential, function of the classical 
field, is defined as the Legendre transform of W 

U eff m = ^W(j)-j<S>, (38) 

where the source has to be seen as a functional of the classical field, by inverting the relation 
( 1361) . The first derivative of the effective potential is 

dU^ ±dWd L _d L 
rf$ V dj d$ d$ J J ' { } 

and the second derivative is then 

dHJ^ l = _d L = _(d ± \- 1 = _J^ >0 m 
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and is positive, for any value of the field $. Therefore, the effective potential of the theory 
must be convex: even if we start with a spontaneous symmetry potential with a concave part, 
when dressing the system, quantum fluctuations rub out the classical spontaneous symmetry 
breaking feature and lead to a convex effective potential. 



3.2 Equivalence with the Wilsonian effective potential 

(I would like to thank Janos Polonyi for the argument given here.) 

The effective potential defined previously, as the Legendre transform of W, is the non-derivative 
part of the proper graphs generator functional, and is denoted U prop . We show here that it is 
equivalent to the Wilsonian effective potential, denoted Uwu s - We note that for this argument 
to be valid, it is essential that we work in Minkowski space time, so as to be able to express 
the Dirac distribution in eq. (|42|) in terms of its Fourier transform. 
For a constant IR configuration $, the Wilsonian effective potential is defined by 

exp {iVU WiU ($)) = J V[<f>\ exp {iS[Q + 0]) , (41) 

where S is the bare action defined at some cut off A, and the dynamical variable which 
is integrated out has non-vanishing Fourier components for \p\ < A. One can also write the 
previous definition as 

exp(zWW$)) = Jv[4>]exp(iS[<f>})8(Jj-V<^ (42) 
= J V[<f>) exp (iS[(f)}) J exp jij Qf <f> - V$) } , 

where J x denotes the integration over space time, j is a real variable, and Jj denotes the 
integration over j. This leads to 



exp (iVU Wils (<$>)) = j ]z(j)exp(-i jj<$> 
= j ]exp(iW(j) -i jj$ 



Jex P (WU prop ($)). (43) 



In the last expression, the integration over j leads to a multiplicative constant, as $ is fixed. 
Disregarding the ^-independent terms, we then obtain 

U W U.($) = Uprapi®), (44) 

which shows the equivalence between the two effective potentials. 



3.3 Towards a flat potential 

Since both definitions of the effective potential are equivalent, Wilsonian transformations to- 
wards the IR, which describe how a system gets dressed by quantum fluctuations, should build 
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a convex potential from a symmetry breaking bare potential. It is indeed the case, and we 
explain here the mechanism. 

The earliest version of exact renormalization group equations was derived in [13J, where a sharp 
cut off is used in order to eliminate degrees of freedom. In this work, an effective theory de- 
scribed by the action Sk~sk, for degrees of freedom with momenta \p\ < k — 5k, is obtained 
from a theory described by Sk, for degrees of freedom with momenta \p\ < k, by integrating 
out degrees of freedom with momenta k — 5k < \p\ < k. An exact evolution equation for the 
coarse grained action Sk with the energy scale k is then obtain in the limit 5k — > 0. In the 
situation of a symmetry breaking bare potential, integrating out degrees of freedom becomes 
problematic when fluctuations with typical momentum k, and above the background 0, satisfy 

k 2 + C/£(0) = 0, (45) 

where is the coarse grained potential defined at the scale k, and a prime denotes a deriva- 
tive with respect to 0. This is the spinodal instability, corresponding to the cancellation of 
restoration force in the system: the propagator of fluctuations with momentum k and above 
the background cannot be defined. Classically, the inverse propagator k 2 + U'^ are {(f)) can then 
be negative when k decreases, leading to an ill-defined path integral. 

An extension of this approach has shown how the Wilsonian flow copes with this spinodal in- 
stability [5J, where it was realized that, when the spinodal instability occurs, non-trivial saddles 
points appear in the functional integral defining the elimination of degrees of freedom. The 
consequence is that k 2 + U'^(4>) actually remains frozen to 0, for any value of k down to the IR 
limit k — 0. The effective potential U e ff = t4 =0 is then flat since Uj! =0 (<p) = 0. This mecha- 
nism corresponds then to a "tree- level renormalization" [14] : it is similar to the famous Maxwell 
construction in Thermodynamics, and is also consistent with other works using a Wilsonian 
approach [15]. An important point is that the flat effective potential interpolates the minima 
of the bare potential, such that the spinodal region goes beyond the range of values of for 
which UH are (4>) < 0. More precisely, if we start with the bare potential defined at some cut off 

UM-^-p + fr, (46) 

where m 2 > 0, it was found numerically in [6J that the coarse grained potential which matches 
the bare potential outside the spinodal region is quadratic: 

f4(0) = ^(2m 2 - k 2 ) - h -<p 2 if 2 <% 2 -A; 2 ), (47) 
2# 2 g 

and its evolution with k is sketched in figUl with the replacement 

k 2 -> t{R - y 2 ) (48) 

Note that, for the sake of clarity, we neglect loop corrections to the potential outside the 
spinodal region. 

To conclude, as the energy scale k decreases, the potential describing the system at this scale 
k gradually rubs out the symmetry breaking features, and in the IR limit becomes completely 
flat between the minima of the UV potential. 
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